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ABSTRACT 



We investigate the type III radio bursts and X-ray signatures of accelerated electrons in a well observed solar flare in order to find the 
spatial properties of the acceleration region. Combining simultaneous RHESSI hard X-ray flare data and radio data from Phoenix-2 
and the Nanjay radioheliograph, the outward transport of flare accelerated electrons is analysed. The observations show that the start- 
ing frequencies of type III bursts are anti-correlated with the HXR spectral index of solar flare accelerated electrons. We demonstrate 
both analytically and numerically that the type 111 burst starting location is dependent upon the accelerated electron spectral index and 
the spatial acceleration region size, but weakly dependent on the density of energetic electrons for relatively intense electron beams. 
Using this relationship and the observed anti-coiTelation, we estimate the height and vertical extent of the acceleration region, giving 
values of around 50 Mm and 10 Mm, respectively. The infeiTed acceleration height and size suggest that electrons are accelerated 
well above the soft X-ray loop-top, which could be consistent with the electron acceleration between 40 Mm and 60 Mm above the 
flaring loop. 

Key words. Sun: flares — Sun: radio radiation — Sun: X-rays, gamma rays — Sun: particle emission 



1. Introduction 

Accelerated electron beams are believed to be responsible for 
both hard X-ray (HXR) and coherent radio emission during solar 
flares. Upwards travelling electron beams propagate through the 
coronal plasma and sometimes escape into interplanetary space. 
Emission from such beams is often observed as coronal and 
interplanetary type III radio bursts. Electron beams travelling 
downwards with small pitch angles enter the dense plasma of the 
chromosphere and are generally seen through bremsstrahlung 
emission in HXR. Before entering the chromosphere, down- 
wards propagating electron beams may also produce reverse type 
III bursts. Despite this wealth of electromagnetic beam emission 
from accelerated electrons propagating in plasma, the location 
of the electron acceleration site and its spatial characteristics are 
poorly known. 

Indirect evidence of electron acceleration sites first came 
from broad band radio spectral obs ervations of pairs of type 
III and re verse type III bursts (e.g. lAschwanden et al] Il995al: 
lAschwan den & Benz 1997). The starting frequencies of these 
burst pairs are found between 220 - 910 MHz, implying a range 
of electron densities in the acceleration region between 6x10^- 
10"^ cm""* for fundamental emission or 1 .5 x 10** - 3 x 10^ cm"-' 
for harmonic emission. These densities are lower than ones ob- 
served in bright soft X- ray loops (2 x 10"' - 2 x 10" cm"^) 
suggesting that the acceleration region lies above the soft X- 
ray loops, being located for example in a cusp reconnection site. 
HXR observations have also been independently used to provide 
insight into typical electron acceleration region heights above 
the photosphere. Through electron time-of- flight analysis us- 
ing H XR emission in the range 20-200 keV (lAschwanden et all 
Il998h . height estimates have been found in the range 20-50 Mm. 



The spatial size of the acceleration region still remains largely 
unknown. 

The simultaneous observation of HXR and metric/decimetric 
radio emission is commonplace during flares and the relation- 
ship between type III bursts and hard X-ray emissions h as 



been studied for many years (see for example iPick & Vilmer 
120081 for a review). The first studies performed by iKane 
(Il972h found a good similarity between HXR and type III ra- 
dio emission, suggesting the two emissions are produced by 
electrons originating from a common acceleration site. Many 
subsequent studies have specifically dealt with the associa- 
tion of coherent type III radio emiss i on and HXR bursts , 
bofli statisticaUy (e.g. [Kane 197 l 119811: iHamilton et alj|1990l: 
lAschwan den et al. 1995a; Arzner & Benz 2005) and for Individ 
ual e vents (e.g. , Kane et al. 1 982; Benz et al . 1983; Dennis et al 
1981 
200C: 



Raoult et al 



Vilmer et al. 



19851: lAschwanden et a l.lll995bl:lRauhn et al 



2002 ^. A more recent statistical study of 201 



flares above GOES class C5 (Benz et al. 2005) reports an 83% 
association rate with coherent radio emission, within the range 
between 4 GHz and 100 MHz. These results suggest that practi- 
cally all flares with HXR GOES class > C5 are associated with 
some form of coherent radio emission. 

An in depth statistical stu dy was carri ed out between radio 
type III bursts and HXRs by ( lKanelll98ih . The study reported 
that the X-ray/type III correlation increases systematically with 
the intensity of HXR and radio emission, the peak spectral hard- 
ness of HXR em ission and the type III burst starting frequency. 
[Hamilton et al.l (fT990) similarly reported the systematic increase 
of HXR/type III correlation with increasing emission intensity 
and to a lesser extent with spectral index of HXR emission. 
Hamilton et al. (1990) also reported a statistical correlation be- 
tween the peak HXR and type III intensities. To produce a harder 
(smaller spectral index) HXR photon spectrum requires a harder 
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electron beam spectrum. A hard electron beam spectrum is an 
attractive attribute for type III producing electron beams as it 
makes it easier and faster for the bump-in-tail instability to oc- 
cur. Faster instability onset ties in very well with the HXR/type 
III correlation increasing for bursts with a higher starting fre- 
quency. 

A temporal correlation bet ween HXR and ra dio pulses has 
been found statistically ( Aschwanden et al.ll995ah where the av- 
erage time delay betwen the HXR pulse and radio pulse start- 
ing frequency was < 0.1 s. Temporal correlations have also 
been found in ma ny individual event studies. Of these stud- 
ies the results by iDennis et al.l (1198 4) find a temporal cor- 
relation with a similar m agnitude to the statistical study by 
lAschwanden et al.1(ll995ah . This, together with previous correla- 
tions, suggests a common acceleration region with HXR produc- 
ing electron beams having either slightly less distance to travel 
or slightly more energetic electrons responsible for the emission. 
The simultaneous analysis of HXR and type III radio bursts is 
thus an attractive diagnostic of flare associated electron acceler- 
ation. 

A few previous studies have attempted to infer properties 
(both height and size) of the common electron acce leration re- 
gion fr om simultaneous HXR and radio observations. iKane et al.l 
used an inferred spectral index from HXR emission to 
estimate the minimum distance required for the type III produc- 
ing electron beam to become unstable. With an assumed elec- 
tron acceleration height of 20 Mm, an altitude of 100 Mm was 
deduced for the starting frequency, in good agreement with the 
spectral observations. Unfortunately the type III frequencies in 
this analysis were too high with respect to the Nangay radio- 
heliograph frequencies for radio imaging at this time so it was 
not possible to c onfirm the starting height of the radio emis- 
sion. iBenz et aTl d 19831) also considered the minimum distance 
required for a type III producing electron beam to become un- 
stable. By modelling both the HXR and radio producing electron 
beam as a Maxwellian they found a weak correlation between 
the type III starting frequency and the electron temperature de- 
rived from HXR observations above 26 keV. However, in the 
event considered a change in electron temperature cannot fully 
account for the initial change in type III starting frequency. The 
authors thus conclude that a movement of the acceleration site 
occurred for this event. 

In another studv' Aschw anden et al.l (1199 5a') uses the assump- 
tion of a common acceleration region producing upward and 
downward electron beams to estimate acceleration times and in- 
fer acceleration region sizes. No frequency gap was observed 
between type III and reverse type III emission so their starting 
frequency separation distance was constrained by the detector 
resolution. Using an assumed density model, this distance was 
found to be < 2 Mm. The instability distance for the electron 
beam was then equated to twice th is resolution. By using a sim- 
ilar analysis of iKane et al.l (Il982l) with observed HXR spectral 
index, acceleration times are found with Af < 0.3 - 3 ms. The 
size of the acceleration region is then inferred at 0.7 Mm. This 
constraint of acceleration site size and times is heavily depen- 
dent upon the assumption of bidirectional electron beams start- 
ing with a separation unresolved by the detector and thus should 
be treated with care. 

Non-thermal HXR emission detected above 20 keV pro- 
vides information on the overall spectrum of accelerated elec- 
trons, usually represented by a power-law in velocity and/or en- 
ergy space. It also provides insight on the density of the non- 
thermal electrons required to explain the HXR emission. Radio 
type III bursts are believed to be formed via coherent plasma 



emission m echanisms requiring the presence of Langmuir wave 
turbulence (iGinzburg & Zhelezniakovll 19581) . The high level of 
Langmuir waves infers a local electron beam-plasma instabil- 
ity. Assuming that some acceleration mechanism in the corona 
produces a power-law in energy space the required unstable dis- 
tribution can be formed due to electron transport parallel to the 
magnetic field. Therefore, the criteria for beam instability and 
hence starting frequency of the type III emission are linked to the 
properties of the acceleration region; the spatial size, the spectral 
index of energised electrons, and the height in the corona. 

In some events upward electron beams are able to travel from 
the corona to the inner heliosphere and create interplaneta ry type 
III bu rsts. Previous work ( Kontai- & R eid 20Q9; Re id & Kontaj 
i2010h has modelled electron beam transport from the Sun to 
the Earth taking int o account the self-consistent generation 
of Langmuir waves. iKontar & ReidI (2009) found that wave- 
particle interactions are required to explain both electron beam 
onset times a nd spectral properties at 1 AU. Developing this 
model further, ^eid & Kontaii(l2010h found the inclusion of den- 
sity turbulence suppressed the growth of Langmuir waves which 
led to a broken power-law spectra similar to observational results 
(fKrucker et al.ll2009l) . 

Electron beams escaping into the inner heliosphere can also 
be detected in-situ near 1 AU and their numbers ha ve been corre- 
lated to the number of HXR producing electrons ('K rucker et al.l 
2007). A correlation is found between the spectral indices of 
both electron populations as well as between the numbers of 
HXR producing electrons and escaping electrons for prompt 
electron events. This again suggests that the X-ray producing 
electrons and the electrons moving upward in the corona origi- 
nate from a common acceleration site. Furthermore, it is found 
that the number of escaping electrons is on average only 0.2% of 
the HXR-producing electrons above 50 keV. 

In this paper we show that simultaneous radio and X-ray ob- 
servations can be used with a simple analytical model to diag- 
nose not only the location but the size of the acceleration region 
as well as the location. This provides one of the first observa- 
tional estimates of acceleration region size. The present analysis 
of the April 15th 2002 at 08:51 UT suggests the acceleration 
region is located at heights of x 50 Mm and has a spatial size 
of a; 10 Mm. Numerical simulations are used with these inputs 
to help validate the results and explore unknown electron beam 
parameters. 

2. Basic concepts: starting frequency of type III 
bursts 

The aim of the following theory is to relate known observational 
variables from flares to unknown flare parameters. Specifically 
we will relate the starting frequency of type III bursts and the 
spectral index of the inducing electron beam to the height and 
size of a flare acceleration region. Initially let us consider the 
propagation of a flare accelerated electron cloud with starting 
characteristic size d located at an initial height r - 0. The initial 
distribution is given by 

fiv,r,t^O)^go(v)expi-\r\/d) (1) 

where goiv) is the initial electron velocity distribution. Langmuir 
waves will be generated when the quasilinear growth rate 
yiv, r) of electrons resonantly interacting ((jjne( r) - kv) 
with Langmuir waves (Dmmmond & Pines 1962; V edenov et al.l 
1962) becomes stronger than the colhsional absorption jc by the 
Maxwellian plasma 
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y(v, r) = — — V — > jc. 



Jc = — 71— 



(2) 



where In A is the Coulomb logarithm, taken near 20 for the pa- 
rameters in the corona, ajp^ir), tieif) and vje are the background 
plasma frequency, density and thermal velocity respectively. 

The initial velocity distribution of solar flare accelerated par- 
ticles is usually taken as a power-law goiy) ~ v^". This distribu- 
tion at f = is stable and will not lead to generation of Langmuir 
waves. At later times f > the propagation of particles leads to 
the formation of a positive slope of the distribution function in 
velocity space. The distribution function changes in time due to 
the propagation (in the case of no energy losses) 



f(v, r, t) = go(v)exp(-\r - vt\/d) 



and the growth rate for Langmuir waves becomes 



7T0J,Jr) 2 

r(v) = — — V 



(3) 



(4) 



We can observe strong Langmuir wave growth remotely via ob- 
servations of type III radio emission. Langmuir wave growth 
should occur at the distance Ar = htypein - hacc from the original 
location where hiypeiii is a height corresponding to the starting 
frequency of type III bursts and hacc is the acceleration region 
height. The distance Ar can be found by equating the growth 
rate for Langmuir waves given in Equation (|4} with y^. giving 



Ar = d\a+ —(vgo(v)) 



(5) 



The quantity vg(){v) - rtt, where «/, is the inducing electron beam 
density. By assuming a coronal background electron density of 
lO' cm"^, a coronal electron temperature of 2 MK and a beam 
density of 10"* cm"-' the second term jcneif) / 7Ta>pe{r)nh is around 
10"^ a. Thus the relation between the acceleration site prop- 
erties and the type III starting frequency is determined mostly by 
the spectral slope of the electron beam giving the simple relation 



^typelll 



— da + ha 



(6) 



The unknown parameters hacc and d can be found from the 
known parameters a and h,ypcin- A key advantage of this rela- 
tion lies with its lack of dependence on the poorly known elec- 
tron beam number density necessary to produce ty pe III emis- 
sion. T he method is similar to what was discussed in lKane et al.l 
(Il982h . The key diff'erence is that we consider an instantaneous 
inject ion at f = with a spatially broad distribution function 
while iKane et alJ ( fT98l considers a temporal injection from a 
point source. 



3. Observations and data analysis 

3.1. Instruments 



In this study, RHESSI <Lin et alJl2002h data is used to analyse 
the HXR emission from downward propagating electron beams. 
By assuming a common acceleration site for upward and down- 
ward electron beams, the X-ray emission observed by RHESSI 
can serve to provide the spectral index for the above predic- 
tion. RHESSI can also provide spatial information regarding the 
HXR. 



When observing the radio emission in the corona, the spec- 
tral information was obtained by using Phoenix-2 radiospec- 
trometer (Messmer et al. 1999) within the frequency range 100 - 
700 MHz. This provides information we use to find the starting 
frequency of the type III radio emission. The spatial information 
for type III radio e mission was found by u sing the Nangay radio- 
heliograph (NRH) (Ke rdraon & DelouiilT997l) in the frequency 
range of 164 - 432 MHz. The spatial radio information observed 
by the NRH can make sure the observed radio flux is emitted 
from the same active region in the solar atmosphere responsible 
for the HXR emission. 



3.2. Selection 

We aimed our study at cases observed simultaneously with the 
NRH, Phoenix-2 and RHESSI. We started our selection from 
a list of events presented in previous observ ationa l analysis 
between coherent radio and HXR emission s dArzner & Benzl 
120051: iGrigis & Ben"3l2004t iBenz et al.ll200l . Of the 58 events 
considered, 10 were found to have coherent radio emission in 
the frequencies covered by the NRH observations. We selected 
one event which had a simple spatial configuration at all NRH 
radio wavelengths, clearly defined starting frequencies and a 
strong HXR flux. This event was chosen as an illustration of the 
method. 



3.3. Measurements 

The spatial overview of the April 15th flare is presented in 
Figure[T]using RHESSI, NRH and SOHO/EIT. The X-ray source 
was imaged using RHESSI in the energy range between 15 and 
30 keV. The higher energies had too few photons to make a rea- 
sonable image above the background noise. The radio source 
was imaged using the NRH in frequency bands from 164 to 
432 MHz and the size increases with decreasing frequency. This 
can provide an estimate regarding the radial magnetic field ex- 
pansion locally in the corona. However, the decrease of spatial 
resolution with decreasing frequency using the NRH has to be 
considered. This decrease behaves as D/A where D is the max- 
imum antenna spatial separation and A is the wavelength of the 
radio emission. The SOHO/EIT 195 image displays informa- 
tion about the overlying plasma configuration, conferring insight 
into the magnetic field structure where the flare originates. The 
temporal evolution of the flare is presented in Figure |2] using 
Phoenix-2, NRH and RHESSI data. 

The spectral index of the X-ray emission, y, was obtained 
using RHESSI spectral analysis of the photon flux /(e) ~ e 
The photon spectral index y was estimated using a power-law fit 
every 2 seconds (half-rotation of the spacecraft). The one sigma 
error associated with the power-law fit was used as the spectral 
index error estimates. 

The starting frequency of the type III radio emission was de- 
termined from the Phoenix-2 data. We used Phoenix-2 data with 
a 1 sec temporal resolution. The mean value of the radio flux 
on the quiet 3 minute interval 08:56 UT to 08:59 UT was used 
to quantify the background level for each frequency channel. A 
threshold of twice this background level was then used at every 
moment in time to determine the starting frequency of the ra- 
dio emission. We then averaged the starting frequency over the 2 
second RHESSI interval. The mean width of the radio channels 
between 100 and 700 MHz was 9.2 MHz so we took 10 MHz as 
the one sigma error on the starting frequency. 
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Fig. 1. The morphology of the April 15, 2002 solar flare. 
Background is SOHO/EIT 195 image. The small red contour 
lines at the base of the plasma loops on the left correspond to 
HXR photons imaged by RHESSI in the 15-30 keV range. The 
large contours on the right hand side correspond to NRH ra- 
dio images at frequencies 432 MHz (blue), 327 MHz (orange), 
236 MHz (pink), 164 MHz (yellow). 




08:52 08:53 
UT Time 



The combined determination of starting frequencies and 
spectral indices was done on two time periods between 08:51:20 
^ 08:51:36 UT and 08:52:38 ^ 08:53:00 UT Both periods 
have a HXR non-thermal spectral index below 7.5 at all points 
in time. Moreover, throughout both periods there is significant 
radio emission above the threshold frequency. 

The temporal evolution of j and the type III starting fre- 
quency is overplotted on the Phoenix-2 data in Figure |3] for the 
time periods defined above. The photon spectral index displays 
an anti-correlation with the type III starting frequency. A clearer 
visualisation of this anti-correlation can be seen when both ob- 
servables are plotted against each other (Figure 0. They have a 
Pearson correlation coefficient of -0.65 due to the starting fre- 
quency decreasing as the spectral index increases. This correla- 
tion suggests that the two variables are related by a linear fit. 

3.4. Acceleration region 

To infer the characteristics of the coronal acceleration region 
from Equation (|6]l we must use some assumptions to obtain 
htypein and a from the type III starting frequency and y. To relate 
the starting height of the type III emission h,ypeni to the starting 
freq uency w e have used the exponential density model derived 
in iPaesold et al.l (l2001i) which assumes second harmonic emis- 
sion for a reference height of around 1 .5 Rs for 160 MHz emis- 
sion. To obtain the electron beam spectral index in velocity space 
from the photon spectral index in energy space, the thick target 
model (Brown 1971; Brown et al. 2006) was assumed. The elec- 
tron beam spectral index a can then be calculated from the pho- 
ton spectral index y through the simple relationship a = 2(y+ 1). 



Fig. 2. Time evolution of the radio and HXR fluxes for the April 
15, 2002 solar flare between 08:51 and 08:54 UT. The top panel 
is the Phoenix-2 radiospectrometer data on a log scale between 
the frequencies 160 and 700 MHz. The middle panel is the 
Nangay radioheliograph flux time profiles observed at 5 discrete 
frequencies from 164 to 432 MHz. The bottom panel are the 
RHESSI HXR counts/second at the three energy ranges 6-12, 
12-25 and 25-50 keV. 

The efi'ect of phot osphere albedo is ignore d as the flare is located 
close to the limb (Kontar & JefiTey"2010|)- 

A positive correlation between the electron beam velocity 
spectral index and the starting height is observed with a Pearson 
correlation coefficient of 0.62 (Figure|5]l. To investigate the cor- 
relation and obtain estimates of the acceleration region prop- 
erties a linear fit w as applied to the data. The routine mpfitexy 
(lMarkwardtll2009h was used to obtain a fit to the data includ- 
ing observational error (Figure |5)- Using Equation (|6) the lin- 
ear fit infers the acceleration region height and size values of 
hacc - 52 + 21 Mm and d - 10.5 ± 1.6 Mm respectively. The 
larger percentage error of the height in relation to the size can 
be observed in the extremes of the fit shown in Figure |5] If this 
linear relationship is statistically significant the slope has to be 
greater than zero. We can test the null hypothesis that the slope 
is zero using the Students t-statistic (e.g. Press et al. 1992). A 
t-score of 6.56 is found with 19 degrees of freedom. Using a 
confidence level of 0.01 we can comfortably reject the null hy- 
pothesis and say the linear relation is statistically significant. 

The radio threshold frequency used to constrain the start- 
ing frequencies had a minor effect on the results if changed 
within reasonable parameters. Different levels (1.5 - 2.5x back- 
ground level) changed the acceleration region properties by 
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Fig. 3. HXR spectral index and frequency spectra of the type III 
burst for two different time periods in the April 15, 2002 event. 
The starting frequencies are plotted as red stars connected by 
dashed lines. The HXR spectral indices are plotted as 2 second 
green bars with error bars in the middle of their integration time. 
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Fig. 5. Scatter plot of electron beam spectral index in velocity 
space, a, calculated from the HXR photon spectral index vs. the 
distance above the photosphere associated with the starting fre- 
quencies of the type III burst. The one sigma observational eiTors 
on both spectral index and height are shown. The green dashed 
line is a linear fit to the data including observational error with 
the green dotted lines showing the extremes of the fit. 

4. Numerical simulations of Langmuir wave 
generation 

To explore these predictions for acceleration height and size we 
can use numerical simulations of electron beams and induced 
Langmuir waves in the solar corona. The simulations allow us to 
validate the observational result given the known initial condi- 
tions and check the validity of Equation Moreover, it allows 
us to explore some of the unknown parameters such as beam 
density and the level of Langmuir waves required for radio emis- 
sion. 




3 4 5 6 7 

HXR Photon Spectral Index 

Fig. 4. Scatter plot of the HXR photon spectral index vs. the 
starting frequencies of the type III burst. The one sigma obser- 
vational errors on both spectral index and starting frequency are 
shown. 



around 5 - 10 %. Higher threshold frequencies caused higher hacc 
and lower d with the converse being true. Threshold frequencies 
< 1.5 or > 2.5x background level caused unrealistic acceler- 
ation region parameters as starting frequencies were either not 
detected or were detected at high frequencies not corresponding 
to the visually observed type III emission. 



4. 1 . Electron beam transport 

The evolution of accelerated electrons can be considered self- 
consistently using weak turbulence theory where we have also 
taken into account binary collisions of energetic electrons with 
the surrounding plasma. We consider the time evolution of an 
electron distribution function /(v, r, f) and the induced Langmuir 
wave spectral energy density W{v, r, t). The energetic electrons 
travel in a beam-like stracture along the magnetic field which 
radially expands outwards into the corona, allowing a one di- 
mensional treatment of their propagation via 



d 



dt (r + ro)2 dr 



(r + roff 



4;rV Q T^gf 
rrP- dv V dv 



lnA|-4 (7) 



dW 
~dt 



d(jJL dW dupeir) dW nojpeir) ^ df 

— y ^ 

dk dr dr dk n^ir) dv 



-ijc + JlW + e^ojpAr)vf\n—. (8) 

VTe 



The background plasma is assumed to be a Maxwellian dis- 
tribution with thermal velocity vje, density nc(r) and plasma 
frequency Upeir). The background plasma damps through col- 
lisions both the electron beam (last term of Equation dTjl) 
and the Langmuir waves (y^. = 7rnp(r)e"^ln A/(m^v^^)) (e.g. 
Melrose 1980). Coulomb collisions also spontan eously induce 
Langmuir waves (last term of Equation dS)) (e.g. iHannah et al.l 
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l2009l) . The first terms on the right hand side of Equations 
(I7I8I I describe the resonant interaction w„f.(r) - kv of electrons 
and L angmuir waves (Drummond & Pines 1962; Vedenov et al. 

This resonant interaction is also responsible for Landau 
damping of Langmuir waves with phase velocities near vje. 



described through yi = 
(iLifshitz & Pitaevskiilll98lh 

As the density of the background plasma changes with dis- 
tance from the Sun, Langmuir waves experience wave refraction. 
This effect is described through the third term in Equation dS) 
where the level of background plasma inhomog eneity is char- 
acterised thr ough / - LOpe{rY^dupe{r)ldr (e.g. lRvutovlll969t: 
lKontai]l2Q01h . 

The spatial propagation of electron beam and Langmuir 
waves are described by the second term in Equations (I7]8) . 
Electron beam spatial propagation also includes the ma gnetic 
field expansion within the corona dReid & Kontajl2010l) . This 
expansion in the corona is modelled through an expanding cone 
which has an angle 9 and starts at length ro from the acceler- 
ation region. The coronal expansion conserves the total num- 
ber of electrons such that for scatter-free propagation, + 
r(xi^ne{r)dr - const along the magnetic field lines. 

4.2. Initial conditions 

The system of equations is solved numerically for an initially 
localised electron beam 



f(v,r,t ^0) ^ go(v)exp\-^-^\, go(v') = ^^^^ 

a / v„, 



(^)" . (9) 



The initial distribution is characterised by d the acceleration re- 
gion size, \',ni„ the minimum electron velocity considered, a the 
beam spectral index, and ni, the density of the electron beam 
above Vmm- The background plasma absorbs Langmuir waves 
with phase velocity close to the thermal velocity via Landau 
damping therefore we set v„„„ to Svjf and safely ignore veloci- 
ties lower than v„„„. 

The t hermal spectral energy density of Langmuir waves (e.g. 
iHannah et al. 2009 ) is the initial condition given by 



Wniv, r,t = Q)^ — — log 

4k^ \ Vje 



(10) 



where is the background plasma temperature set at 2 MK and 
Kb is Boltzmann constant. This is the expression for the thermal 
level of a Maxwellian plasma when collisions are weak. 

The background dens ity Wg(r) was defined using the same ex- 
ponential density model dPaesold et al.ll200l]) that was assumed 
in the previous section. 

4.2.1. Observational constraints 

The values derived from the observations in the previous sec- 
tion constrain some of the key input parameters for the simu- 
lations. The starting height, hacc - 52 Mm which corresponds 
to a background density of »g - 3 x 10^ ^"^^ us ing the expo- 
nential density model given in lPaesold et al.l (1200 ih . This gives a 
plasma frequency of 500 MHz relating to second harmonic emis- 
sion of 1000 MHz. The characteristic beam size d - 10.5 Mm. 
The HXR spectral index y is found from the RHESSI observa- 
tions (Figure|5]l which allows us to constrain the electron veloc- 
ity spectral index as 6 < a < 16. 



The NRH images of the type III radio emission allow us 
to observe how the radio source increases with decreasing fre- 
quency. Such an observation can provide information regarding 
the magnetic field expansion. The size of the radio emission at 
237 MHz is approximately twice the size of the radio emission 
at 432 MHz taken at 30% of emission level. This was mea- 
sured around the two peak times of emission at 08:51:21 UT 
and 08:52:42 UT. However, the wavelength A is approximately 
twice as large at 237 MHz compared to 432 MHz and so the 
angular resolution of the NRH is increased by two. Moreover, 
scattering by density inhomogeneities will increases the appar- 
ent size of the coronal radio source more at higher wavelengths 
(lBastianlll994 ). We thus in the present case cannot observa- 
tionally resolve any significant radial expansion of the magnetic 
field. Such a scenario is equivalent to type III producing electron 
beams propagating along thin coronal structures as observed in 
lTrottetetal.1 (e.g. il982l) : iPick et all (e.g. 120091) . or having very 
small radial expansion of the magnetic field in the low corona. 
The expansion is much smaller than what would be expected 
for the inner heliosphere, where the magnetic field expands as a 
cone with an angle around 40° (e.g. Steinberg et al. 1985). 

The density of electron beams responsible for type I II emis- 
sion is believed to be small with .Krucker et all (l2007l) finding 
them 0.2% of the density of the downward propagating elec- 
tron beams responsible for HXR emission above 50 keV. With 
an initial background density of 3 x lO** cm""* providing the up- 
per limit to the downward propagating electron beam, the up- 
ward propagating electron beam was injected with a density of 

- lO'* cm"-' above 11 keV. However, observations show time 
dependent intensities of HXR photons which is related to the 
density of the inducing beam. Such results could indicate the po- 
tential need to consider a changing beam density. We note that 
Equation (|6]l is independent of the density of the electron beam. 
The starting frequencies found from the upwardly propagating 
electron beam should thus be insensitive to rather large changes 
in beam density. 

4.3. Numerical Results 

A high level of Langmuir waves is required to induce type III 
emission. We can estimate the starting height, /i/v;,^///, from the 
simulations through the ratio of Langmuir wave spectral energy 
density to its initial thermal level W{v,r,t)/WTh{v,r,t = 0) or 
W/Wth- The first point in phase space when W/Wth exceeds a 
certain level can give us insight into how electron beams with 
different spectral indices become unstable. 

The numerical results are presented in Figure |6]for a variety 
of different W/Wth levels. By assuming Langmuir waves pro- 
duce radio emission when they reach a certain level of W/Wth, 
we can treat the curves in Figure |6] in a similar manner to the 
observational results. By applying a linear fit to each curve, we 
can obtain an estimate of the initial simulated acceleration region 
height and size using Equation As we know the actual initial 
simulated values for ha^c and d, these estimates allow us to check 
how accurate the method is for obtaining good estimates. Such 
a fit also provides a numerical check for the simplified analyti- 
cal relation Equation (|6]l represents. We find the closest fit to the 
simulated hacc and d comes from the line where W/Wth = 10^ 
giving hacc =43.5 + 5 Mm and c/ = 12.4 ± 0.6 Mm. These vari- 
ables are within 20% and 15% of the original numerical values 
respectively. 

The results in Figure |6] show a small variation between the 
heights corresponding to 10^ < W/Wth < 10^. Provided there 
are enough electrons to generate sufficient Langmuir waves for 
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Fig. 6. Heights corresponding to high levels of Langmuir waves 
from an unstable beam with density lO'* cm"-'. Symbols and 
colours correspond to different levels of Langmuir wave growth. 
Low levels (10 W/Wrh) correspond to spontaneous emission of 
waves. High levels correspond to beam-plasma instability. 



radio emission, a change in the beam density has minimal effect 
on the starting height hnpeiii- Increasing or decreasing the beam 
density by one order of magnitude changed the inducing height 
of W/Wti, - 10^ by at most 14 % when a - 6 with a mean 
over all spectral indices of 3 %. Changing the beam density will 
only vary the level of Langmuir waves which are produced upon 
the electron beam becoming unstable. This result confirms the 
density independence of Equation (|6]l where instability of the 
electron beam is mainly dependent upon the spectral index and 
size of the electron cloud. 

The ratio W{v, r,t)/WThiv, r,t - 0) also provides informa- 
tion regarding the Langmuir wave phase velocities and onset 
times when the waves exceeds certain thresholds. As electrons 
resonantly interact with Langmuir waves, the phase velocity of 
the Langmuir waves conveys information regarding velocities of 
the inducing electrons. The phase velocities corresponding to the 
points in Figure |6] get smaller as the electron spectral index in- 
creases (softer spectrum). An example of the velocity variation 
is presented in Figure |7]for the level W/Wti, = 10^. Similarly the 
time required for the beam to induce Langmuir waves at a cer- 
tain level increases for larger spectral index (Figure 0. Beams 
with larger spectral indices (with the same beam density) have 
less high energy electrons. We thus expect Langmuir wave emis- 
sion to be induced by lower energy electrons which take longer 
to become unstable. 

We can explore how the inclusion of an expanding magnetic 
field would change the simulation results. By ignoring the in- 
crease in radio source size at lower frequencies from NRH reso- 
lution and scattering we can assume the doubling of size between 
432 and 237 MHz is because of a radially expanding magnetic 
field. Using the assumed density model the magnetic field would 
then expand as a cone with an angle of = 6°. By then assum- 
ing an acceleration site size which is as wide as it is long, we can 
constrain ro in EquationQto be 30 Mm below the solar surface. 
Running the simulations and using the same method described 
above gives estimations of hacc and d which are 12% and 40% 
of their original values. The expanding magnetic field causes the 
electron beam to induce a high level of Langmuir waves further 
away from the acceleration site. The frequencies corresponding 
to these heights do not agree as well with the observed starting 
frequencies of the type III bursts. 



2.0 



1.5 



H 1.0 



a 
O 



0.5 



0.0 



8 10 12 14 

Velocity Spectral Index 



16 



18 



Fig. 7. Top: electron velocity spectral index plotted against the 
first phase velocity at which Langmuir waves exceed the thresh- 
old W/Wtii = 10^. Bottom: electron velocity spectral index plot- 
ted against the onset time required for the electron beam to in- 
duce Langmuir waves exceeding the threshold W/Wrh - 10^. 



5. Discussion and conclusions 

We used simultaneous observations of radio and hard X-ray 
emission during a solar flare to gain insight about the accelera- 
tion region of energetic electrons. With a simple model we have 
shown through an analytical relation how the starting height of 
type III emission and the spectral index of the electron beam 
can be related to the height and vertical extent of an accelera- 
tion region. By combining HXR spectral information with the 
starting frequencies of the type III bursts, we have derived for 
our event an estimate on the acceleration site height and size of 
hacc = 52 ± 21 Mm and d = 10.5 + 1 .6 Mm respectively. We have 
also used self-consistent numerical simulations of an electron 
beam which can induce Langmuir waves in a background coro- 
nal plasma. The simulations checked our predicted acceleration 
region values and allowed us to explore unknown parameters of 
the electron beam and Langmuir wave distributions. 

The value found for hgcc agrees with values in the range 
20 - 50 Mm, deduced from electron time-of-flight analysis for 
HXR emission (iAschwandenll2002h . This scenario indicates an 
acceleration region in the corona well above where SXR are 
imaged. The error on hacc is quite large but within the 95 % 
range of 2-sigma the acceleration region remains within the 
corona. The value found for d is roughly a n order of magni- 
tude hi gher than previously found before in lAschwanden et al] 
(II995ah . Assuming the relation in Equation ^ the acceleration 
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size in lAschwanden et al] (Il995ah would not be able to pro- 
duce significantly varying starting frequency of type III emis- 
sion given a static acceleration site . Such a small acceleration 
site predicted by Aschwanden et al.l (Il995a ) may be relevant for 
type III radio bursts when very little or no evolution of the start- 
ing frequency can be observed. 

Using the estimates for h^cc and d we ran self-consistent 
numerical simulations of an electron beam able to resonantly 
induce Langmuir waves in the background coronal plasma. 
We analysed the distance required for a large magnitude of 
Langmuir waves to be induced through beam instability for a 
variety of initial beam spectral indices. A linear fit to the ini- 
tial beam spectral index and the height associated with strong 
Langmuir wave production gave a good estimation of the ini- 
tial acceleration region height (within 20%) and the initial ac- 
celeration region size (within 15%) . The result fits with the an- 
alytical predictions from Equation (|6]l and hence the relation is 
a powerful diagnostic tool for flare acceleration site properties. 
In Une with the analytical equation the electron beam instabil- 
ity criteria was significantly dependent upon a and d and almost 
independent on the beam density, which was confirmed by nu- 
merical simulations. The simulations also gave an estimate of 
W/Wth ^ 10^ as the magnitude of Langmuir waves required to 
produce coherent radio emission. The discrepancies found in the 
acceleration region properties are due to additional terms present 
in the numerical simulations which were not present for the sim- 
ple analytical expression. 

It is also possible to explore how a different assumed ini- 
tial electron beam distribution in space will affect our results. 
Initially in Equation [1] we assumed a tent distribution for the 
electrons in space. We now consider an initial electron beam dis- 
tribution which is Gaussian distributed in space such that 

f(v,r,t^O)^go(v)exp(-^/d\ (11) 

Th e instability criteria f or this distribution was already discussed 
in iMel'Nik & Konta3 (Il999h and, assuming small collisional 
damping, gives the relation 

^2d^f^ + hacc■ (12) 

The dependence of htyp^ju on the square root of the spectral 
index originates from the in the exponential for the elec- 
tron distribution function. Unfortunately the observational er- 
rors on electron beam spectral index were too large to distin- 
guish between the two models (Eq ([TJ and Eq (fTTT i). Even with- 
out observational error estimates such a fit to the data gives 
hacc = 33 + 51 Mm and d - 22 + 1 Mm which is not defined. 
More detailed observations are thus needed to differentiate be- 
tween different initial electron distributions in space. 

Another assumption we considered was a static acceleration 
site during the entire event. iKane & Raoultl (119811) considered 
a moving acceleration site which decreased in altitude to ex- 
plain why a type III burst's starting frequency increased with 
time. Assuming the magnetic nature of reconnection, any source 
movement would typically be at the Alfven veloci ty. At heights 
around 100 Mm, this is typically around 1 Mm s ' ( Arregui et al.l 
l2007h . The Alfven velocity is too slow to account for the vary- 
ing starting frequencies of the type III emission observed in the 
April 15th flare considered. Moreover, the acceleration region 
would have to move upwards and downwards to account for the 
evolution of the starting frequencies. Our results do not rule out 
the acceleration region moving in altitude but this will probably 
not be the dominant process for determining dynamic type III 
starting frequencies on a time scale of seconds. 



Flares associated with the same active region responsible for 
the 08:51 UT flare on the April 15th 2002 have been analysed 
previously. Sui & Holman (2003) found a coronal HXR source 
above the loop-top HXR source during another flare around 
23:00 UT on the same day. The high coronal HXR source was 
initially detected at an altitude of 25 Mm and moved with a 
speed of 0.3 Mm s ' up to an altitude of 40 Mm as the HXR 
flux increases. Moreover, the higher energy photons (15-16 keV) 
are detected at lower altitudes than the low energy photons 
(10-11 keV). This is indicative of an electron beam streaming 
down from a high acceleration region with high energy elec- 
trons having a larger stopping distance than low energy elec- 
trons. Such a scenario fits with the derived high acceleration 
region hacc ~ 50 Mm we found in this study. A similar result 
was found for other high coronal sources (' Liuetal.ll2008ll200l 
where high energy photons are imaged at lower altitudes than 
low energy photons. 

In conclusion, we stress that simultaneous HXR and radio 
observations are a tool to estimate the otherwise unmeasurable 
sizes of the acceleration site. The results from our first trial of 
the relation given by Equation |6] suggest that this size can be 
a: 10 Mm located at height ^ 50 Mm, occupying a substantial 
fraction of the corona. The size is larger than the HXR sources 
which are typically observed with RHES SI and in the rang e be- 
tween a few Mm up to a few tens of Mm ( E mslie et al .'2003). We 
hope to continue this avenue of research by extending it to more 
events to build a statistical picture of flare acceleration charac- 
teristics. Future studies should also have a higher flux of HXR 
to better constrain the deduced electron beam spectral index. 
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